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The zero-temperature phase diagram of p-orhital two-component fermionic system loaded into a one¬ 
dimensional optical lattice is mapped out by means of analytical and numerical techniques. It is shown that 
the p-band model away from half-filling hosts various competing superconducting phases for attractive and re¬ 
pulsive interactions. At quarter filling, we analyze the possible formation of incompressible Mott phases and 
in particular for repulsive interactions, we find the occurrence of a Mott transition with the formation of fully 
gapped bond-ordering waves. 

PACS numbers: 71.I0.Pm, 75.I0.Pq 


I. INTRODUCTION 


Ladder systems have been the focus of much theoretical and 
experimental work over more than two decades. One theoret¬ 
ical motivation was to investigate the dimensional crossover 
between the well-known one-dimensional (ID) physics and 
the two-dimensional (2D) case, in the search of 2D non-Fermi 
liquid physics. A second reason stems from experiments and 
the study of ladder compounds, such as the famous telephone 
number one SiT 4 _a;Caa;Cu 2404 i which has a superconducting 
phase at high pression and for a small hole density. ' 

The simplest ladder model is that of a two-leg ladder, made 
of two coupled fermionic chains. In stark contrast to the sin¬ 
gle chain case, the two-leg ladder system displays a super¬ 
conducting phase with d-wave superconductivity for repulsive 
interactions which stems from the doping of a spin-gapped 
Mott insulating phase at half-filling.^’^This gave the belief that 
the two-leg ladder problem already contains seeds of the rich 
physics of the cuprates. 

Two-leg ladders have thus become over the years a funda¬ 
mental system for the study of low-dimensional strongly cor¬ 
related fermions. Various exotic quantum phases have been 
predicted theoretically depending on the form of the coupling 
between the two chains and the filling."^^^*' Experimental re¬ 
alizations of two-leg ladder systems are clearly called for, in 
particular, to investigate the rich physics of the weak-coupling 
regime. In this respect, ultracold fermionic gases are a promis¬ 
ing way to study two-leg ladder problems thanks to the high 
level of control on interchain hopping and interactions. 

The ladder geometry might be created by considering double¬ 
well optical lattices for instance.^^~^^ 

A second possible way is to load a two-component Fermi 
gas in a optical lattice and consider higher-lattice orbitals, 
typically the p-band, to simulate a fermionic two-leg ladder 
system. More precisely, we consider, in this paper, a two- 
component Fermi gas which is loaded in a ID optical lat¬ 
tice (running along the z-direction) with moderate strength of 
(harmonic) confining potential muP‘{x^ -h p^)/2 in the direc¬ 
tion (i.e., xy) perpendicular to the chain.^®“^*^ It is assumed 


that all the s-level of the oscillator are fully occupied while 
the p-level, i.e., Px,y, are partially filled. The resulting lattice 
fermionic model has been derived in Refs. 26-28 within the 
tight-binding approximation and takes the following form: 

.^p-band — f ^ ^ i^rna, i-\-l T 
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where m = px,Py is the orbital index and a = t)i is the 
’’spin” index or internal components of the underlying cold 
atoms. In Eq. (1), m = Y.ma describes the den¬ 

sity operator at site i and a pseudo-spin operator for the orbital 
degrees of freedom has been defined: 



na, z? 


( 2 ) 


where a = x,y,z are the Pauli matrices and a summation 
over repeated indices is implied. Model (1) can be viewed as 
a two-chain fermionic system without interchain hopping but 
with density-density interchain interactions and pair-hopping 
processes between the two orbitals (see Fig. 1). 

The continuous symmetry of model (1) for general Ui ,2 is: 
U(l)c X SU(2)s, where U(l)c denotes the U(l) (charge) sym¬ 
metry related to the conservation of the total number of atoms, 
while SU(2)s is the internal global SU(2) (spin) symmetry of 
the two-component Fermi gas. We note that the U(l)o contin¬ 
uous symmetry for the orbital degrees of freedom is explicitly 
broken in Eq. (1). However, the use of harmonic potential 
in the xy direction implies a constraint of the two coupling 
constants [/i 2 and the investigation of model (1) along the 
harmonic line Ui = 3 U 2 enjoys an U(l)o symmetry corre¬ 
sponding to rotation along the y-axis in the orbital subspace. 
The U(l)o symmetry is trivially realized when U 2 = 0 where 
the model reduces to two decoupled SU(2) Hubbard chains. 
Finally, the Ui = U 2 case displays also an U(l)o symmetry 
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(1) is then deduced in the weakly-interacting regime for in¬ 
commensurate filling and at quarter filling with Fermi vector 
kp = (oq being the lattice spacing). In Sec. Ill, DMRG 
calculations are carried out to investigate the intermediate and 
large-coupling regimes to fully determine the phase diagram 
of the p-band model. We present then our concluding remarks 
in Sec. IV. 


FIG. 1. (color online) The two-leg fermionic ladder representation 
of the p-band model (1). The orbital degrees of freedom plays the 
role of the legs of the ladder. Besides the kinetic term with am¬ 
plitude t for each orbital, there are local interactions at each site i: 
density-density intra-orbital (amplitude (7i, in red) and inter-orbital 
(amplitude U 2 , in green), orbital exchange (amplitude U 2 , in blue) 
and inter-orbital pair-hopping that breaks U(l)o generically (ampli¬ 
tude U 2 , in orange). 


since it is directly related to t /2 = 0 after a redefinition of the 
orbital pseudo-spin operator. We consider, in this paper, the 
most general case where (7i 2 are not fine-tuned for two main 
reasons. On the one hand, departure from the harmonic line 
Ui = 3 U 2 can be investigated by breaking the axial symme¬ 
try of the 2D harmonic trap with the introduction of quartic 
potentials.^* On the other hand, the harmonic line has been 
wrongly identified in Refs. 26 and 27 and we want to make 
contact with their results obtained when U 2 = 4(7i/9. 

At half-filling, it has been shown that model (1) along 
the harmonic line with repulsive interaction describes a Hal¬ 
dane phase^® of the spin-1 Heisenberg chain.The latter 
is known to be the paradigmatic ID topological phase pro¬ 
tected by symmetry.*® For an attractive interaction, a Hal¬ 
dane phase for the charge degrees of freedom*'’** has also 
been predicted in Ref. 27 by means of a strong-coupling ap¬ 
proach and density-matrix renormalization group (DMRG) 
calculations.** In Ref. 28, it has been shown that there is adi¬ 
abatic continuity between weak and strong coupling regime 
and these two Haldane phases are related by a spin-charge in¬ 
terchange symmetry. The half-filled p-band model (1) along 
the harmonic line thus paves the way to realize experimentally 
non-trivial Haldane phases in the context of cold fermions. 

In this paper, by means of a low-energy approach and 
DMRG calculations, we map out the zero-temperature phase 
diagram of the p-band model (1) for incommensurate filling 
and at quarter filling which best avoids three-body looses. In 
stark contrast to the conclusion of Ref. 27 for attractive in¬ 
teraction, we show that the p-band model (1) away from half¬ 
filling does not behave as an attractive ID Hubbard chain with 
the formation of a 2kp CDW but instead is a representative 
of the physics of two-leg fermionic ladders. In particular, we 
find that the hallmark of its phase diagram is the emergence of 
various competing superconducting phases for attractive and 
repulsive interactions as well as a Mott transition with the for¬ 
mation of fully gapped bond-ordering waves (BOW). 

The rest of the paper is organized as follows. In Sec. II, 
we present our low-energy approach based on a mapping onto 
Majorana fermions and a one-loop renormalization group 
(RG) analysis. A zero-temperature phase diagram for model 


II. THE LOW-ENERGY APPROACH 

In this section, we perform a weak-coupling approach to 
model (1) away from half-filling. We focus on incommensu¬ 
rate filling and a special commensurate filling where kp = 
2 ^. The starting point of this approach is the introduction of 
left-right moving Dirac fermions from the continuum limit of 
the non-interacting Hamiltonian (1) with C7i 2 = 0: 

Cla,i (3) 

with X = iao, I = Px,Py, and a =t,In stark contrast to the 
half-filled case of Ref. 28, there is a ’’spin-charge” separation 
which strongly simplifies the analysis and the full Hamilto¬ 
nian of the p-band model (1) decomposes into two commuting 
pieces; 


n=Ha + 'Hs, (4) 

with \Hc^T~Ls\ = 0, where Tfc governs the physical proper¬ 
ties of the charge degrees of freedom and TTs refers for all 
remaining non-Abelian degrees of freedom, i.e., spin, orbital 
and spin-orbital degrees of freedom. 


A. Charge degrees of freedom 


Let us first consider the simplest part in Eq. (4), i.e., 
the charge degrees of freedom. For incommensurate fill¬ 
ing, on general grounds, the charge part takes the form of a 
Tomonaga-Luttinger Hamiltonian;*’* 




Vc 

2 




(5) 


where Vc and Kc are the Luttinger parameters. In this low- 
energy approach, the charge excitations are described by the 
bosonic field 4)c and its dual field ©c. The explicit form of the 
Luttinger parameters for the p-band model can be extracted 
from the continuum limit (3) after standard calculations. We 
find; 


a /1 + gd^vp 

Vc = Vp^/l+ gc/TTVp, (6) 

where vp = 2tao sin(kpao) is the Fermi velocity and pc = 
2 ao(Ui + 112 ). 

For incommensurate filling, no umklapp term appears and 
the charge degrees of freedom display metallic properties in 
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the Luttinger liquid universality class.^’^ However, at quarter 
filling, umklapp processes might be generated leading to an 
additional term in Eq. (5). Such perturbation can be found 
by means of higher orders of the RG calculations or from 
lattice symmetries.In this respect, the one-step translation 
symmetry expresses as follows in terms of the charge bosonic 
field; Tag : $c ^"0 + \/ 7 r /4 . The charge perturbation 
with the smallest scaling dimension which is compatible with 
the translation symmetry is then: Vc = —gu cos(-\/l67r$c). 
The full charge Hamiltonian density for the commensurate 
filling kp = 4 ^ becomes then equivalent to the quantum 
sine-Gordon model at = IfiTr; 


He = 


Kr 


- gu cos(a/1(^T>c). 


(7) 

This model is exactly solvable and the development of the 
strong-coupling regime with a charge gap Ac cotTesponds to 
Kc < 1 / 2 , while for Kf. > 1/2 the charge degrees of freedom 
remain gapless in the Luttinger universality class. Using the 
estimate ( 6 ), we find the position of the critical line for a Mott 
transition within the bosonization approach: 


Ulc + U 2 C — 
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( 8 ) 


such that Kc = 1/2. One has to be very careful with this 
estimate since the full expression of Kc as a function of C 7 i .2 
( 6 ) is strictly valid only in the weak-coupling regime. One 
needs a complementary approach, i.e., numerical calculations 
to conclude on the existence or not of a Mott transition for the 
lattice model (1). Note that the mechanism of the Mott tran¬ 
sition is very similar to the one which occurs in the repulsive 
U(4) Hubbard chain^^’^® or in the problem of spin-3/2 cold 
fermions. 

For Kc < 1/2, the charge degrees of freedom are fully 
gapped and the strong-coupling regime of the sine-Gordon 
model (7) is described by the pinning of the charged field on 
the minima: 

($c) = > 0 

(®c) = \A74, gu < 0, (9) 

p being an integer. The sign of the umklapp coupling constant 
gu is difficult to fix but we expect 9u ^ 0 from an higher- 
order perturbative expansion in the weak-coupling regime so 
that:($c) = The DMRG calculations of Sec. Ill will 

shed light on the position of the pinning of the charge field 
from the determination of the Mott-insulating phases of the p- 
band model. The low-lying excitations are massive kinks and 
antikinks which interpolate between the ground states of the 
quantum sine-Gordon model (7). The charges associated to 
these charge excitations are 


Q = zL 



dx dx^c = ±1, 


( 10 ) 


in units of the electron charge. We have thus massive holon as 
low-lying excitations in the charge sector. 


B. Non-Abelian sector 

The low-energy physics of the remaining degrees of free¬ 
dom, included in Hs of Eq. (4), can be inferred from a map¬ 
ping onto six Majorana fermions as in Refs. 14-16, 41, and 
42 

A simple way to obtain this correspondence is through 
the introduction of four chiral bosonic fields = 

Px,Py',<T = t; i from the bosonization of Dirac fermions:^’^ 



where the bosonic fields satisfy the following commutation 
relation; 


[^laRj^l'a'L] = (12) 

The presence of the Klein factors ensures the cotTect an¬ 
ticommutation of the fermionic operators. The Klein fac¬ 
tors satisfy the anticommutation rule {kict, Ki'a'} = ‘^Su/Saa’ 
and they are constrained so that T^ = 1, with T = 

Hereafter, we will work within the T = 1 

sector. 

The next step of the approach is to introduce a bosonic ba¬ 
sis which singles out the different degrees of freedom, i.e., 
charge, spin, orbital, and spin-orbital degrees of freedom:"^* 

%^^L,R = 

%^iL,R = 

%,^L,R = ^ ^ 

<^PyiL,R = ^ (4>c - - 4>o + (13) 

From these new bosonic fields, one can now consider a 
refermionization procedure by introducing six left and right 
moving Majorana fermions through: 

^R + = -^= exp {iV^<S>sR^ 

= -J= exp - /V4//4>ol) 

exp 

^R + exp (/V47r4>soi?), (14) 

where 771 , 2 , 3,4 are again Klein factors which ensure the ade¬ 
quate anticommutation rules for the fermions. Using this cor¬ 
respondence, the interacting part of Hs can be expressed in 
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terms of these Majorana fermions: 



+ + ^G^Ril^R^L + ^7^R^t^R^L^ ( 15 ) 


with the identification: 


Ai 

A 2 

A 3 

A 4 

As 

As 

Ay 


—ao{Ui 

+ U 2 ) 

—2aoU2 


o-o{U2 — 

Ui) 

0 


2 aoU2 


0-0 {Ui — 

U 2 ) 


0 . 


Using standard calculations, the one-loop RG equations for 
model (15) can de derived: 

A2 = —A1A2 + 7;—A3A7 + 7—A4A6 

TT ZTT ZtT 

A3 = —A1A3 + 7—A2A7 + 7—A4A5 

TT ZTT ZTT 

A4 = —A1A4 + 7—A2A6 + 7—A3A5 

TT ZTT ZTT 

3 1 

As = ^A3A4 + —AsAy 
3 1 

As = —A 2 A 4 + —AsAy 

3 1 

Ay = 77^2A3 -f ^AsAg. ( 16 ) 

As often in ID, the RG equations enjoy some hidden dis¬ 
crete symmetries: 


fit : A2,3,4 —y —A2,3,4 

AI 2 ■ A2,6,y —y —A2,6,y 

^^3 : A4^5_6 —)■ —A4_s,6 

^4 : A3_5_y —> —Ag^s.yi (17) 

which correspond to the existence of duality symmetries on 
the Majorana fermions"^^ 

fit : el’"’' ^ 

Ci 

^ 3 : ci 

: Ci ^ , (18) 

while the right-moving Majorana fermions remain invariant. 
These chiral transformations for each Majorana fermion cor¬ 
respond, in the continuum limit, to the well-known Kramers- 
Wannier duality symmetry of the underlying ID Ising model 
in a transverse field. The latter transformation maps the or¬ 
dered phase to the disorder one which is ordered in terms of 


the disorder operator.^ Starting from the ordered phase of the 
six underlying Ising models in a transverse field, the transfor¬ 
mations (18) lead then to four possible different gapful phases. 
A numerical RG analysis of Eqs. (16) is necessary to deter¬ 
mine which phases is indeed reached starting from the initial 
conditions (16). 
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FIG. 2. (color online) Top: Numerical phase diagram obtained by the 
one-loop RG analysis of Eq. (16). Bottom: For each phase denoted 
I-IV, we give the dominant superconducting correlations (see text) as 
well as the duality symmetries which relate the phases between them. 
We display the central charge c = 2 or c = 4 which characterizes 
the different quantum phase transitions of the model. 

We did this numerical analysis and, as depicted in Fig. (2), 
it reveals that the RG flow goes in the strong-coupling regime 
in the far infrared (IR) along four asymptotic lines: 

/ : A2 = A3 = A4 = —A5 = —Ag = —Ay = —Ai 
II: X2 = —A3 = —A4 = —Ag = Ag = Ay = —Ai 
III : A2 = A3 = —A4 = —Ag = —Ag = Ay = Ai 
IV : A2 = —A3 = A4 = — Ag = Ag = —Ay = Ai. ( 19 ) 
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Along these rays, the interacting part of the effective Hamil¬ 
tonian (15) simplifies as follows; 




ei 






i I > 


t 07^6 


( 20 ) 


with p > 0. Using the duality symmetries (18), all these mod¬ 
els reduces to the same model which takes the form of the 
SO(6) Gross-Neveu (GN) model with interaction;"^^ 

^gn = | ■ (21) 

This is an example of a dynamical symmetry enlargement 
(DSE) where an SO(6) symmetry, i.e. a higher symmetry than 
the initial symmetry of the p-band lattice model, emerges at 
low-energy. This phenomenon occurs in a large variety of 
models with marginal interactions in the scaling limiU"^ as, for 
instance, in the half-filled two-leg Hubbard model'® and in the 
SU(4) Hubbard chain at half-filling,"'^ where an SO(8) sym¬ 
metry occurs at low energy. The emergence of an SO(6) sym¬ 
metry has also been obtained within the one-loop RG analysis 
of various models of doped two-leg ladders. 


C. Phases of the p-band model 

The main interest of this SO(6) DSE scenario stems from 
the fact that the isotropic RG ray is described by a massive 
integrable field theory (21) when g The develop¬ 

ment of the strong-coupling regime in the SO(6) model leads 
to generation of a non-perturbative fermionic mass, i.e. the 
emergence of a spin-gap phase where all spin, orbital and 
spin-orbital excitations are fully gapped. The nature of the 
underlying electronic phases of the p-band model can then 
be inferred by a straightforward semiclassical approach of the 
SO(6) model and the application of the duality symmetries 
(18). Alternatively, one can perform a direct bosonic semi¬ 
classical approach of the different models in Eq. (20) using 
the identification (14). In the following, we proceed to this 
analysis by identifying the phases of Eig. 2. 


1. Phase I: Hi duality 


that the RG flow is attracted in the far IR along a special ray 
(I) of Eq. (19). This phase is described by the Hi duality and 
the resulting physical properties of that phase are governed by 
the interacting Hamiltonian of Eq. (20). A straightfor¬ 
ward semiclassical approach of the latter model leads to the 
following pinning of the bosonic fields; 


($,)=0, = (©,<,) =0. (22) 

The electronic properties of phase I, which includes the har¬ 
monic line Ui = 3 U 2 , depend also on the charge degrees of 
freedom that are decoupled. 

Let us first consider the Kc >1/2 case. One has then a gap¬ 
less charge excitation whereas all remaining spin, orbital, and 
spin-orbital excitations are fully gapped. All 2kp densities are 
short-ranged due to the pinning of the dual spin-orbital field 
©so- Interestingly enough, the leading electronic instability 
in this phase turns out to be a bond-centered superconducting 
instability with order parameter; 


^BCSbond — d + y {jPx Py)^ (23) 


which is odd under the Z 2 orbital exchange symmetry (px O 
Py). A similar superconducting instability has been intro¬ 
duced in the study of doped two-leg electronic ladder in 
Ref 48. In the continuum limit, we have; 


^BCSbond — Sin( ^ “f ^Pyi^Px't (24) 
~ ^Pxi^Pyt ~ ^Pyt^Pxi) ■ 

One can then bosonize this operator by means of the identifi¬ 
cation (11) and the physical basis (13); 

4sin(A:iBao) 


O, 


BCSbond — 




— iy/^Qc 


(25) 


x[ — i sin s/n^s cos s/tt^o sin s/nQs 
+ cos s/n^s sin y/n^o cos y/n<ds 


Taking account of the vacuum expectation values (22), we get 
(^BCSbond — so that the equal-time correlation function 

of the pairing operator has a power-law decay; 


(C^lcsbobd (^) acsbbbd (0)) ^ (26) 


When 1/2 < Kf. < 1, the bond-centered BCS superconduct¬ 
ing instability is strongly enhanced with respect to the non¬ 
interacting case. We have thus a dominant superconducting 
instability for repulsive interaction. 

While the 2kp densities are short-ranged in phase I, Akp 
densities might compete with the superconducting instability 
(23). Many Akp density terms can be written in the continuum 
description like ~ L]^RioL]_^Ri-a for instance. 
The latter gives in the bosonized language; 

Piil ~e®^''=cosv^<i>o, (27) 

so that in the spin-gapped phase (22), one has ~ 

gi%/47 r$c^ and therefore a power-law decay for its correlation 
function; 

(pIK ( 0 )) 


In the first phase, when (7i > 0 and Ui > t /2 of Eig. (2), 
the numerical analysis of the one-loop RG equations reveals 


(28) 
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Since > 1/2, one observes that the Akp density correla¬ 
tion function (28) decays much faster than the superconduct¬ 
ing one (26). In summary, the leading instability of phase I 
with Kc > 1/2 is the bond-centered BCS superconducting 
one. 

In the Mott-insulating phase with K,. <1/2, all excitations 
are fully gapped. The bond-centered BCS superconducting 
instability (23) is now short-ranged since the charge held 
is pinned (9) and thus its dual held 0^ is a strongly huctuating 
held. In contrast, the Akp density of Eq. (27) can have now a 
non-zero expectation value, breaking spontaneously the trans¬ 
lation symmetry. We thus expect a Mott phase with a two-fold 
degenerate ground state at quarter hlling kp = ^. The na¬ 
ture of this phase depends on the sign of the coupling constant 
gu of the umklapp perturbation in model (7). In this respect, 
we introduce the following order parameters as in the study of 
the extended quarter-hlled two-leg Hubbard ladder:^^ 

C^BOW,r,qy(*) = (~1) + 

cr 

+ COSQy Cp^a-.i+l^PyC^A + 

(T 

+ cos qy (29) 

with qy = 0, TT. One can obtain a continuum and bosonized 
representation for these order parameters. The Akp contribu¬ 
tion of the latter plays a crucial role and can be derived using 
the results of Refs. 34, 49, and 50 and we get for the leading 
contribution: 

(-l)*^c^i+iC/,T,i + H.c. ~ sinV47r($c + 

<7 

(-1)* ^ ~ cos V47r($c + ei^o), (30) 

<7 

with = 1, Cp^ = —1. From these results, we hnd the 
bosonized descriptions of the order parameters (29): 



cos \/Pr^c sin 


OBOW,r,o 

sin •s/Urfl’c cos •\/47r$o 



sin \/47r$c sin V47r$o 


C’cDW,r,o 

cos y/Pt^c cos a/TttT'o. 

(31) 


Taking into account that in phase I, we have ($o) = 12, 

we have either a uniform Akp BOW or uniform Akp CDW 
depending on the sign of the umklapp term in Eq. (7): 

(ObOW„,o) 7^ 0; if ffu < 0 

(OcDW„,o) 7^ 0; if 9u > 0. (32) 


2 . Phase II: 0,2 duality 


In the second phase, when U 2 > 0 and Ui < C/ 2 , the one- 
loop RG flow is now attracted in the far IR by the asymptote 
(II) of Eq. (19). This phase is described by the 02 duality. 
A straightforward semiclassical approach of model Ti.^ in 
Eq. (20) leads to the following pinning of the bosonic fields 
in this phase: 

($,)=0, (0o)=O, ($«o)=0. (33) 

In this phase, as we will see in Sec. Ill, numerical results 
And that Kc > 1/2. All spin, orbital, and spin-orbital ex¬ 
citations are fully gapped in this phase and the charge de¬ 
grees of freedom are gapless since Kc > 1/2. Due to the 
presence of attractive interactions, one may call this gapless 
phase a Luther-Emery phase as in the spin-1/2 Hubbard chain 
with U < OP Its physical nature is very different from that 
of phase I since the expectation values of the bosonic fields 
(33) are different. As it can readily be seen from Eq. (33), 
we have again the condensation of a dual held, here the or¬ 
bital one 0o, which implies that the 2kp-CDW operator is a 
strongly fluctuating order. As in phase I, the only possible 
CDW quasi-long range order in this phase is a Akp CDW. 
However, the CDW of Eq. (27) becomes now short- 
range since the orbital dual held 0o condenses in phase 11. 
In this respect, we have to consider another Akp CDW insta¬ 
bility: Pikp E<t Ll^^jRp.PlyaRpya- The latter can be 
directly expressed in terms of the bosonic held: 

PaL ^ cos . (34) 

Using the expectation values (33), we deduce the leading 
asymptotics of the equal-time correlation function of the Akp 
CDW operator: 

(35) 

As it will be seen in the next section, we have Kc >1/2 and, 
one does not expect, on general ground, that this d/cp-CDW 
will be the dominant instability in this phase. As in phase I, 
a superconducting instability turns out to be the leading one. 
The latter is defined by the following order parameter: 


^BCS— — Cpy^ ,iCpy]^^i, (36) 


which is odd under the Z 2 orbital symmetry (p^. o py) and 
antisymmetric with respect to the spin degrees of freedom 
(to/,). The order parameter (36) can be expressed directly 
in terms of the bosonic fields: 


Gbcs- —-e 


(37) 


cos cos cos 



We expect in the weak-coupling regime that < 0. A uni¬ 
form Akp BOW is thus stabilized which is two-fold degener¬ 
ate and breaks the translation symmetry. The latter phase is 
similar to the BOW phase obtained in the quarter-filled spin- 
3/2 SO(5) chain model.:’® 


Using the expectation values (33), we get Obcs- — 
so that the equal-time correlation function of the pairing oper¬ 
ator (36) has a power-law decay: 

(O^CS- (a:)OBCS- (0)) ~ 


(38) 
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which dominates the Akp CDW ones (35) when Kc > 1/2. 


3. Phase III: fla duality 


The next phase, defined by C /2 < 0 and Ui < C/ 2 , is de¬ 
scribed by the asymptote (III) of the one-loop RG flow. The 
harmonic line of the p-band model with attractive interaction 
belongs to this phase which is described by the duality with 
interacting Hamiltonian of Eq. (20). The bosonic fields 
of the bosonization approach are now pinned to the values: 

($,)=0, (0„) = ^, ($,„)= 0. (39) 

A spin-gap is formed and a gapless c = 1 phase emerges in 
this attractive regime with Ui ^2 < 0 since the umklapp term 
cannot gap out the charge degrees of freedom when > 1. 
In close parallel to the previous cases, one can determine 
the nature of the leading electronic instability of this Luther- 
Emery phase by means of the bosonization approach com¬ 
bined with the pinning (39). The dominant Akp CDW is the 
one of phase II, given by Eq. (34) with the power-law behav¬ 
ior (35). The relevant superconducting instability for phase III 
is defined by: 


Obcs-i- = (40) 


which is even under the Z 2 orbital symmetry (p^: ^ Py) and 
antisymmetric with respect to the spin degrees of freedom 
(t^i)- In terms of the bosonic fields, it reads as follows: 




(41) 


i cos sin s/nQo cos s/nAiso 
+ sin s/nA>s cos s/nQo sin v/tt^’s 


Using the vacuum expectation values (39), we immediately 
get C9bcs-i- — so that the equal-time correlation 

function of this pairing operator has a power-law decay: 


(^) ^^BCS+ (0)) (42) 

Since Kc > 1, we conclude that the superconducting instabil¬ 
ity (40) dominates the dfcp-CDW ordering (35). The Luther- 
Emery phase is thus governed by a superconducting instability 
(40) in stark contrast to the 2kp-CDWphase predicted by the 
DMRG study of Ref. 27. The physics of the p-band model with 
attractive interaction is thus not similar to that of ID attractive 
Hubbard model as emphasized in Ref. 27. 


4. Phase IV: H 4 duality 

The last phase of the p-band model corresponds to the re¬ 
gion where U 2 < 0 and Ui > U 2 . The numerical analysis 
of the RG flow shows that, here, the one-loop RG flow is at¬ 
tracted in the far IR by the special line (IV) of Eq. (19). The 


resulting phase is described by the 174 duality with interact¬ 
ing Hamiltonian of Eq. (20) which leads to the following 
pinning for the bosonic fields of the bosonization approach: 

($,)=0, ($o) = 0, (0«o)=O. (43) 

As before, this pinning leads to the formation of a gapless c = 
1 phase when Kc >1/2 where the charge degrees are the only 
critical modes of the problem. We now consider the standard 
d-wave superconducting instability of the two-leg electronic 
ladder to determine the nature of phase IV: 


Oecsd = Cp^t,iCp^i,i - Cp^i,iCpyt,i, (44) 


which is even under the Z 2 orbital symmetry (p^ Py). The 
bosonized expression of this superconducting instability reads 
as follows: 




Oecsd — —^ °' e 


(45) 


i sin sin s/ttA>o sin s/ttQs 


+ cos s/tt^s cos s/n^o cos y/nQg 


From this expression, we observe that this operator is a fluctu¬ 
ating order, i.e., has short-ranged correlation, in the previous 
phases, while, in phase IV, one has from the pinning (43): 
decsd ~ The d-wave superconducting instability 

(44) becomes dominant in phase IV with the power-law be¬ 
havior: 


((!7j5CSd(^)C^BCSd(0))~x-i/2^> (46) 

From the pinning (43), we find that the 2 fc 4 ;’-CDW operator 
is short-ranged while the Akp CDW of phase II, given by 
Eq. (34), has a power-law behavior in phase IV: 

(pfkl (^) difel (0)) ~ , (47) 

with subleading exponent when Kc >1/2 with respect to the 
superconducting instability (46). 


5. Quantum phase transitions 

From the duality symmetries (18), we can, as well, discuss 
the different quantum phase transitions that occur in the p- 
band model by investigating self-dual manifolds. However, 
Fig. 2 reveals that the transitions belong to special lines of 
the lattice model: Ui — U 2 and U 2 — 0. From Eq. (1), 
the 1/2 = 0 describes two decoupled quarter-filled spin-1/2 
Hubbard chains with Ui coupling constant. When Ui < 0, a 
spin-gap is formed and therefore the phase lEphase III tran¬ 
sition is critical with c = 1 -f 1 = 2 gapless charge modes. 
When Ui > 0, all degrees of freedom are gapless since the 
quarter-filled spin-1/2 Hubbard chains is known not to exhibit 
a Mott transition.^ The phase I/phase IV transition is thus crit¬ 
ical with central charge c = 4. The transition between phase 
I and phase 11 is located along the Ui = U 2 > 0 line of the 
p-band model. From the definition of the p-band model (1), 








8 


one can show that the latter line corresponds to two decou¬ 
pled repulsive quarter-filled spin-1/2 Hubbard chains which 
does not have relevant umklapp process. A c = 4 behavior 
should occur along this quantum phase transition. Finally, the 
last transition between phase III and phase IV belongs to the 
U 2 = Ui < 0 line which takes the form of two decoupled 
attractive quarter-filled spin-1/2 Hubbard chains with a c = 2 
quantum critical behavior. All these results can be derived by 
investigating self-dual manifolds of the RG Eqs. (16). Since 
the transition lines are located on special high symmetry lines, 
one would expect the location of these transitions to remain 
universal beyond the weak coupling regime where this anal¬ 
ysis is valid. This will be indeed confirmed numerically in 
Sec. Ill B. A summary of the phases and quantum phase tran¬ 
sitions, obtained from the low-energy approach, can be found 
in the bottom of Fig. 2. 


III. DMRG CALCULATIONS 
A. Determination of the phase diagram 

We will now determine the phase diagram of the N = 2 
p-band model at quarter filling (i.e. 1 particle per site) using 
numerical simulations with the DMRG algorithm. This will 
allow to go beyond the one-loop RG analysis done in Fig. 2 
and check the analytical predictions for realistic intermediate 
or strong couplings (t/i, U 2 ) (we fix t = 1 as the unit of en¬ 
ergy). Moreover, numerical data are needed to determine the 
numerical value of the Luttinger parameter Kc which allows 
to compute the dominant correlations. 

Typically, we have used open boundary conditions (OBC) 
and lengths L = 64 and L = 128, keeping up to 4000 
states when computing correlations in order to keep a dis¬ 
carded weight below 10“^° in most regions, although simu¬ 
lations where 0 < U 2 < Ui were found to be more difficult 
to converge (discarded weight around 10”^, see discussion 
below). For practical purpose, we have mapped the p-band 
model onto an equivalent (pseudo)spin-l/2 (where the pseudo 
spin corresponds to the orbital) fermionic models on a 2-leg 
ladder, and we have implemented the abelian U(l) symmetry 
corresponding to the conservation of particles spin. 

Our main result is presented in Fig. 3 where we plot the 
phase diagram vs (t/i/f, [/ 2 /f) obtained by computing var¬ 
ious correlation functions on L = 128 systems. We will 
present the numerical data below, but we can already discuss 
the different phases. First of all, the (quantum phase) transi¬ 
tion lines are found to ht U 2 — 0 and Ui = U 2 , which is 
expected by symmetry since they correspond to special lines 
of the model, see Sec. IIC 5. A detailed analysis will be given 
in Sec. Ill B. Second, it is remarkable that phases II, III and IV 
found in the weak-coupling analysis are confirmed to exist in 
a wide range of parameters. Last, we have computed the Lut¬ 
tinger parameter Kc using the dominant correlation, which is 
always of some BCS type, in each of these phases. 

Overall, we observe a rather good agreement between the 
numerical phase diagrams obtained by solving RG equations 
(Lig. 2) or the full microscopic model using DMRG (Lig. 3), 


except for the large Mott phase at this commensurate filling. 
Lor incommensurate filling, we predict that the RG phase di¬ 
agram would be identical to the numerical one. 



FIG. 3. (color online) Numerical phase diagram obtained by DMRG 
simulations with length L = 128. Numerical values correspond to 
the Luttinger parameter Kc in each region. 

We will now present some numerical data that were used to 
compute this numerical phase diagram. We have relied mostly 
on computing superconducting coiTelation functions (in vari¬ 
ous channels) as well as density ones. To avoid spurious ef¬ 
fects due to OBC, we have chosen to compute correlations as 

Coix) = {0^{L/4 + x)0{L/A)) (48) 

so that we can deteiTnine if they decay algebraically or are 
short-ranged. In Ligs. 4-5-6, such correlations are plotted and 
correspond respectively to phases II, III and IV. In all cases, 
we emphasize that dominant correlations are found to be the 
superconducting ones, in different channels (see below). In¬ 
deed, the Luttinger parameter Kc is always larger than 1/2 in 
the critical phases, otherwise a Mott phase is stabilized. Some 
values of Kc are given on the phase diagram in Lig. 3. In 
particular, for phases II, III and IV, we observe an adiabatic 
continuity from weak to strong coupling, while phase I on the 
contrary is much reduced and replaced instead by Mott phase 
at intermediate coupling and polarization at strong coupling. 


1. Phase II 

In this region of the phase diagram, our numerical results, 
shown in Lig. 4 for instance when Ui/t = —6 and 1 / 2 /! = 6, 
are in perfect agreement with the RG predictions and domi¬ 
nant BCS- predictions. However, numerics is needed to com¬ 
pute the Luttinger parameter Kc and determine whether pair¬ 
ing correlations dominate over density ones. Using Lqs. (35)- 
(38), our numerical fits for both BCS_ and CDW correlations 
are compatible with a single Kc >1/2, for instance 0.65 for 
the parameters set chosen on the plot. All the other correla¬ 
tions are short-ranged as expected. The leading instability in 
this phase is therefore the BCS_ superconducting one. 
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FIG. 4. (color online) Absolute values of correlation functions (see 
text) vs distance for {Uxjt = —6, Ui/t = 6) and L = 128 obtained 
with DMRG. This corresponds to our phase II. Fitting the dominant 
BCS- correlation with Eq. (38) leads to Kc = 0.65, which is also 
consistent with the TIcf CDW correlation function (see Eq. (28)). 


2. Phase HI 

The phase III region is particularly interesting since it 
could be achieved experimentally using the harmonic trapping 
scheme (i.e. Ui = ’iU^) with attractive interactions. As seen 
in Fig. 5, we have found that the dominant correlations are of 
the pairing type, in the BCS+ channel. In all this region, they 
can be fitted using Eq. (42) to get Kc > 1 (values are given 
in the phase diagram). In such a case, the dominant density 
correlations are uniform and decay as 1/a;^ as expected. 



FIG. 5. (color online) Absolute values of correlation functions 
(see text) vs distance along the harmonic line in the attractive case 
(Uilt = —6, U 2 /t = —2) for L — 128 obtained with DMRG. This 
corresponds to our phase III. Fitting the dominant BCS+ correlation 
with Eq. (42) leads to Kc = 1.45. 

Note that this dominant superconducting correlation func¬ 
tion was not computed in Ref. 27, where only CDW signal 
was discussed. As a result, the dominant instability in this 


Luther-Emery phase is a superconducting one, in stark con¬ 
trast to the 2 /cf-CDW phase predicted earlier.^^ 


3. Phase IV 

In the region corresponding to phase IV, our numerical re¬ 
sults in Eig. 6 indicate that dominant correlations are of the 
pairing type again, but in a different BCS^ channel, as ex¬ 
pected from the low-energy analysis. In the whole region, 
we can fit these correlations using Eq. (46) to extract Kc, or 
equivalently we could fit the subleading CDW correlations 
with Eq. (47), although it would be less precise. Some val¬ 
ues of Kc are given in Eig. 3 and are always larger than 1/2 so 
that no Mott transition is found. 



EIG. 6. (color online) Absolute values of correlation functions (see 
text) vs distance for (Ci/f = 6, U^/t = —6) and L = 128 obtained 
with DMRG. This corresponds to our phase IV. Fitting the dominant 
BCSd correlation with Eq. (46) leads to Kc ~ 0.65, which is also 
consistent with the Akp CDW correlation function (see Eq. (28)). 


4. Phase I 

The region 0 < Ui < U 2 , which should correspond to 
phase I (see Fig. 2) according to RG solution is more involved 
to analyze. We have plotted on Fig. 7 the finite-size scaling of 
the (tt, 0) BOW, i.e. the kinetic energy difference measured 
in the middle of the chain. Data are given along the harmonic 
line Ui = 3(72- Extrapolations are compatible with a finite, 
albeit small, value for intermediate interactions (for instance 
Ui = 6, U 2 = 2), while it seems to vanish in the weak and 
strong coupling regimes. It is known that the weak coupling 
is difficult numerically since the non-interacting starting point 
correspond to two decoupled fermionic chains with a large to¬ 
tal central charge c = 4. Eor instance, the Mott transition in 
the SU(A^) Hubbard model at filling 1/iV has been discussed 
quite extensively to occur for a finite critical {U /t)c when 
N > 2 based on bosonization and quantum Monte-Carlo re- 
sults^^ as well as DMRG ones,^® while some older DMRG 
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simulations had indicated a vanishing Uc = 0.^' Since the 
charge gap is expected to open in an exponentially way, this is 
clearly a difficulty for any numerical technique. However, this 
regime is perfectly suited for bosonization and weak-coupling 
RG; indeed, since the Luttinger parameter Kc = 1 in the non¬ 
interacting case, umklapp processes are irrelevant so that finite 
interactions are necessary to enter the Mott phase. 

Thus, our interpretation is the following: (i) for small in¬ 
teraction parameters, RG analysis should be valid and phase I 
with dominant BCSbond correlations is expected with Kc > 
1 /2 so that umklapp processes are irrelevant; (ii) for interme¬ 
diate interactions, a Mott phase occurs with (tt, 0) BOW and 
exponentially decaying BCS correlations at large distance; 
(iii) for very large interactions (i.e. Ui = 18 and U 2 = 6), 
we have noticed that the ground-state is ferromagnetically po¬ 
larized (hence degenerate): the polarized ground-state can be 
simply understood as two decoupled spinless fermionic chains 
(one for each orbital), which energy is independent of both Ui 
and U 2 - Intuitively, such a state could be stabilized at large 
(7’s since other competing states will have higher energies. 
Still, our numerics is clear on this point as can be checked by 
comparing ground-state energies in different sectors (varying 
the number of particles per spin), or measuring Kc = 1 from 
charge correlations (data not shown). 



FIG. 7. (color online) Finite-size scaling of the (tt, 0) BOW order 
parameter along the harmonic line Ui = 3172. 


B. Nature of the phase transitions 


energies to get more reliable results (see Fig. 11 of Ref. 40 for 
instance). Our numerical results should be compared to the 
analytical predictions made in Sec. IIC 5. 

In Fig. 8(a), we have plotted typical data in phases II, III 
and IV. By removing oscillations using the bond kinetic en¬ 
ergy as an additional fitting parameter, we can obtain smooth 
functions that perfectly agree with the expected behavior (49) 
with c = 1. 



FIG. 8. (color online)(a) Block entanglement entropies vs size x 
obtained for L = 128 and various {Ui, 1 / 2 ). Oscillations have been 
removed using the bond kinetic energies t{x) and fitting parameter A 
(equal respectively to 2, 1.25 and 2.4 (see text). In all cases, we can 
extract a central charge c ~ 1 as expected, (b) Similar data along the 
quantum phase transitions which can be fitted with c = 2.0 (using 
A = 0.325) in the attractive case and with c = 4.0 (using A = 1.5 
and 1, plus a shift for clarity) in the repulsive case using smaller 
L = 64 system. 


A simple way to obtain information on the phases, or on 
their phase transitions, is through the measurements of block 
von Neumann entanglement entropy which is known to scale 
as:^^ 

5'vN = (c/6) logd(a:|L) -F Cst (49) 

where c is the central charge and d{x\L) = (T/tt) sin(7ra:/L) 
the conformal distance. On a finite system with OBC, due to 
Friedel oscillations, the fitting can be more involved and one 
can use for instance the knowledge of the local kinetic bond 


Now, considering the expected quantum phase transition 
along the lines Ui = U 2 and U 2 = 0, we can clearly iden¬ 
tify different behaviors for attractive vs repulsive interactions 
(as expected since there will be a finite spin gap or not re¬ 
spectively). In Fig. 8(b) for attractive interactions, we can fit 
our entanglement entropies data perfectly with c = 2 using 
L — 128. In the repulsive case, our numerical data are not 
converged on the same size keeping up to m = 8000 states, so 
that we plot instead data obtained with i = 64 and m = 8000. 
In this case, we cannot remove entirely the oscillations in the 
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data, but we can get a very good fit using the expected c = 4 (grant 2015-P0677). The authors would like to thank CNRS 

central charge. Overall, we get an excellent agreement with for financial support (PICS grant), 

the theoretical predictions. 


IV. CONCLUSION 

We have presented a comprehensive study of the most 
general model relevant for one-dimensional p-band two- 
component fermionic cold gases with local interactions only. 
We have concentrated here on incommensurate filling and 
quarter-filling. 

Using a state-of-the-art low-energy approach, supple¬ 
mented with a one-loop RG numerical analysis, we have 
found that generically, the charge sector decouples from the 
spin-orbital one which is gapped. As a consequence, most 
of the phase diagram is occupied by standard Luttinger liq¬ 
uid phases with a single gapless charge mode (hence a central 
charge c = 1). Nevertheless, we have clarified the nature 
of the dominant instability and have found that it is always 
of some BCS superconducting kind, in one of the following 
channels: BCSbond. BCS_, BCS 4 . and BCS^, see Eqs. (23)- 
(36)-(40)-(44) respectively. In particular, an interesting bond- 
centered superconducting instability emerges along the har¬ 
monic line for the repulsive interaction. The nature of the 
phase transitions between these four superconducting phases 
is also elucidated and found to behave with central charges 
c = 2 or c = 4. 

Our numerical simulations do confirm that phases with 
dominant BCS_, BCS+ and BCS^ extend from weak to 
strong coupling and occupy large regions in the phase dia¬ 
gram. In particular, for attractive interactions and harmonic 
trapping, BCS+ correlations are the dominant ones, differ¬ 
ent from the 2fci;’-CDW phase predicted earlier.^^ In the last 
region where the bond-centered superconducting instability 
BCSbond is expected at weak coupling, our DMRG data at 
quarter-filling (i.e., one particle per site) indicate that a Mott 
phase intervene with fully gapped bond-ordering waves at in¬ 
termediate coupling, and spontaneous polarization at strong 
coupling. We have also numerically confirmed the nature of 
all the quantum phase transitions present in this model. 

The p-band two-component fermion mode, studied in this 
paper, is thus an interesting model to explore superconduct¬ 
ing instabilities of doped two-leg fermionic ladder systems as 
well as the occurence of a Mott transition. Given the recent 
progress in realizing SU(A^)-symmetric Fermi gases,we 
hope that some of the phases discussed here could be real¬ 
ized experimentally and probed using local spectroscopy tech¬ 
niques or by measuring short-range correlations. 
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